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ABSTRACT 


A method has been described to estimate the effects of finite chemical 
reaction rates on the one-dimensional expansion of air in hypersonic nozzles. 
The calculations have considered stagnation temperatures from 4000 to 6000°K, 
stagnation pressures from 100 to 1000 atmospheres, and a range of wedge and 


axisymmetric nozzle shapes. Vibrational equilibrium is assumed. 


A relaxation length criterion has been applied to a simplified kinetic 
model of air to determine the approximate location of freezing for flow with 
finite reaction rates in each nozzle configuration. Equilibrium composition 
profiles obtained by machine calculation were used for the calculation of the 
relaxation lengths. The resultant chemically frozen expansions have been 
calculated and are presented in tabular and graphical form. In all cases 
freezing occurs fairly early inthe nozzle. Further, freezing in the nozzle 
is delayed by an increase in stagnation pressure, an increase in stagnation 
temperature, and by the use of gradually expanding nozzles. The effect of 
freezing is to reduce the pressure, velocity and temperature at a particular 
area ratio from the corresponding equilibrium values, and to increase the 
density and Mach number. The change in temperature and Mach number 
may be considerable whereas the density and flow velocity are relatively 


unaffected by freezing. 
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NOTATION 


Dimensional 
A cross sectional area of nozzle 
D_ dissociation energy per mole 
Af, change in standard free energy for Joye reaction 
A; enthalpy of Ath species including energy of formation, cal/mole 
& Planck's constant 


if 
H enthalpy of mixture, cal/gm 


Re recombination rate coefficient, Eee jenoie* sec 
Ke equilibrium constant = Re [Re , mole /em? 


Kop Keg equilibrium constant referred to partial pressures and mole 
fractions respectively 


K Boltzmann constant 
£ characteristic nozzle dimension Lia 
Le characteristic throat dimension, ¢m 
m' mass flow_o'u’A’ 

p’ fluid pressure 

f 


Pr local relaxation length, Equation (35) 


rr local relaxation length based on infinite-rate equilibrium, 
Equation (36) 


R, universal gas constant 
Jo entropy (per unit mass) at reservoir conditions 
¢ 
Je molar entropy of <th species at standard pressure 


T local absolute temperature of gas 


vii 


a kp dissociation rate coefficient, ome /mole sec 


us local gas velocity 

% distance from nozzle throat 

%; concentration of species i, moles/gm 

@& characteristic dissociation temperature D/%o 
' characteristic vibrational temperature hp/K 


molecular weight of undissociated mixture 


characteristic frequency of molecule 
fluid density 


2) 
p* 
FP molecular weight of mixture 
y 
/ 
ws 
Pe: 
Y parameter Z£n5 Be [st 


Dimensionless 
Q tangent of nozzle semi-angle 


Ajg member of correction coefficient matrix in Newton-Raphson 
method 


A A'/As 


C number of linearly independent species 


an 


function used in Newton-Raphson method defined by Equations 
(A-2) and (A-3) 


hh; ti: [Ry To" 


hy fractional correction to the #th independent variable in Newton- 
Raphson iteration 


m mass flow pul 


vili 


fri number of chemical elements present 


M Mach number 
M; chemical formula of 4th species 
A % ft 
% exponent in area relation A= 7+ z 
, f 
pp / fe 


q; mole fraction of / th component in a hypothetical system containing 
components only 


Qg number of gram atoms of #th chemical element present 
r r Tt 
er gly a4 


S number of species 


if Sikes 


aya or 1 Ro 


of ? 
jag Sores 
11 [Ro To’ 
«BE 
Koy 


x, mole fraction of (th species 
4Y; formula vector of ( th species 


jh atoms of &th chemical element per molecule of g th chemical species 


wah a , Equation (32) 


Kronecker delta 


maximum permissible fractional correction to mole fractions 


ey ee 
O/T,” 


ix 


stoichiometric coefficients 


Py 5 


Be Wok 


Subscripts 


a, sf) k, 4 sum, product or matrix indices 


oO denotes condition at reservoir 
€ denotes local equilibrium value 
00 denotes infinite-rate equilibrium value 


* denotes value at throat, A= 1 


f denotes value at freezing 


—_———— CZ Ze Maa a sme eee SES we 


I INTRODUCTION 


The stagnation enthalpies encountered in the flow of gases through 
nozzles at low supersonic velocities are usually sufficiently low so that 
chemical kinetics are unimportant. Such flows may be accurately described 
by simple equations which consider the participitation of only the transla- 
tional and rotational degrees of freedom of the molecules (constant 7 }. 
However, in the high temperature flows occurring in hypersonic nozzles, 
such simplified considerations are no longer applicable because of the 
interchange of energy between the internal and translational degrees of 
freedom and the coupling of the chemistry with the gas dynamics through 
dissociation. Figure 1 shows that at stagnation conditions pertinent to the 


present study, a significant portion of the total energy resides in the internal 
energy modes. 


The gas dynamic and chemical behavior of a high temperature gas 
undergoing expansion in complete equilibrium may be calculated with reason- 
able precision (Sec. II). Complete chemical equilibrium, however, implying 
infinite reaction rates, is a condition which may only be approached, to a 
greater or lesser degree, in the expansion of high temperature real gases. 
Ao a peta Sala of chemical nonequilibrium, in particular nozzle 


flows involving only a single dissociating molecular species, have 


shown that the flow properties may depart considerably from their local 
equilibrium values. The effect of finite chemical reaction rates, as well 


as nozzle geometry, may result in actual freezing of the gas composition 
as some point in the nozzle. 


The purpose of this report is to present quantitative estimates for the 
effect of flow freezing on the departure of test section properties from their 
equilibrium values in the expansion of air through hypersonic wind tunnel 
nozzles. The analysis has considered stagnation temperatures from 4000°K 
to 6000°K and stagnation pressures from 100 to 1000 atmospheres, repre- 


senting a range of stagnation enthalpies of either current or future interest. 


The expansions are considered for two nozzle configurations; a two-dimensional 


wedge-type nozzle with a sharp throat and a hyperbolic axisymmecric nozzle 


which is conical for large area ratios. The solutions cover a comprehensive 


range of geometries for both types of nozzle. 


The outline of the report is as follows. The equilibrium chemical 
composition of the gas at the reservoir condition is first computed for all 
stagnation states by the method outlined in Sec. I]. The equilibrium (infinite - 
rate) expansions are then calculated, the solutions providing both the gas 
dynamic behavior and gas composition as a function of temperature. Because 
of their usefulness, the results of the equilibrium calculations are included 


in this report in tabular and graphical form. 


The concept of a local chemical relaxation length” is then introduced 
(Sec. I11) and used to develop an approximate criterion for freezing. This 
concept was first used to develop a freezing criterion in Ref. 3, for both a 
pure diatomic gas or diatomic gas plus any diluent. Some air calculations 
were also performed in Ref, 3. The present kinetic model for air assumes 
finite reaction rates for oxygen, with nitrogen and nitric oxide assumed to 
remain in chemical (infinite-rate) equilibrium. A relaxation length is 
calculated, based upon the equilibrium flow properties, and used to determine 
the approximate nozzle location at which appreciable freezing has occurred, 


i.e., to predict, approximately, the final chemical state of the expanding flow 
for finite reaction rates. 


Finally (Sec. 1V), calculation of the nozzle expansions are completed 
with all chemical species fixed at their ''freezing'' mass fractions. The 
equilibrium solutions serve as a basis of comparison for the frozen (zero- 
rate) solutions. Both the equilibrium and frozen expansion processes are 
isentropic. /* The comparisons have been referred to as providing quanti- 
tative ''estimates" since although the calculations of the frozen expansions 
are in themselves exact, the criterion for the selection of the equilibrium 


state to serve a8 input to the frozen calculations is approximate. 


II SOLUTION OF THE EQUILIBRIUM FLOW 


A. Governing Equations 


It is assumed that the gas is a mixture of ideal gases, that the 
flow is quasi-one-dimensional, and that diffusion, heat conduction and 
viscous processes are of negligible significance. Steady, adiabatic, 


critical flow is assumed to exist in the nozzle. 


The equations which determine the flow are, 


Conservation of Momentum: 
d £ 
udu’ + SE 2 O (1) 


= pressure 


p 
A = density 


Global Conservation of Energy: 


VES od r 
a t+ H =H, (2) 
where H’= mixture enthalpy (cal/gm) 


Global Conservation of Mass: 


¢ 


m= pula (3) 


where Mm mass flow 


r f 
A = cross-sectional area 


The state equations are, 


Thermodynamic State: 


Ro 
B 


ip es (4) 


where & = molecular weight of mixture and is given by 
(5) 
i y 
where J = mass concentration of species ¢ in moles per unit mass. 
Caloric State: 
, = Lf 
H' = a, % A; (6) 


where h; = molar enthalpy of 4th species including thermal and formation 


contributions. h; is only a function of temperature. 


In dimensionless form, the equations are written 


udu pe <0 (7) 
es 

ee ne (8) 
ms: nuh (9) 
p ae (10) 


Zot (11) 


whe re 


: PRS Ts ‘ 
HH 
= p/ po’ 
P* ALP 
ae oy eae 
A= AVA, 
A = molecular weight of undissociated (cold) gas 
R, = Universal gas constant 


Ay = throat area 


Except for the molecular weight, M , the reference conditions for the non- 


dimensional forms are those at the equilibrium reservoir state. 


B. Equilib rium Analysis 


The methods and notation of Brinkley® were used for the equilib- 
rium analysis wherever applicable. At the end of this section, the Brinkley 


technique is extended to include the computation of states along an isentropic 
path. 


The chemical formula of the 4th species may be expressed in 
vectorial form as 


a CG, pgs 208 Kins Mm ) (é #7, 2, ree g) ee 
12 


where &; % is the number of atoms of the #th chemical element per molecule 
of the ( th species and#m is the number of elements present. If the rank of 
the matrix of %;, is c, there exist C linearly independent y; , and(‘-c ) 
dependent Yi which are given by linear combinations of the C independent 
ones, 


Ey yy (ln ert ee 2,08) (13) 


As expressed by Eq. (13), Py are coefficients in the linear dependencies of 
vectors Yo upon a ° 


In the equilibrium analysis the independent species are designated 
( 4 = 1,2.....¢c ) and are called "components.'' The number of components 
C usually equals the number of chemical elements present ™. Equations (13) 
correspond to (f-C) chemical reactions postulated to form the dependent species 


if the formula vectors are replaced with chemical symbols. 


It is useful to consider a hypothetical system. This system has the 
same atomic constitution as the actual system, but is chemically combined 
such that only components are present. If Qp is the number of gram-atoms 
of the &th element present, and Qj is the number of moles of the 7 th compon- 
ent contained in the hypothetical system, then 


Z xy, i = Che 7, Ba a) (14) 


The 95 are normalized by 


ee (15) 
4y £y 


ed 


thereby scaling the total size of the system if $j are considered to be numbers 
of moles. 


Global mass conservation may be expressed by 


¢ 
EL Pa PA z EN eo A z Gi oy (10) 
s aCe 2 


where A; is the molecular weight and 7y is the number of moles of / th species. 


Mass is conserved in each formation reaction of Eqs. (13) so that 


> My LaCie Geese eens s.) (17) 


Substituting “~~ from Eq. (17) and equating coefficients of Fj in Eq. (16) gives 


s 


hj as h; Ne * 9; Gat2, 222) (18) 


Division of Eq. (18) by the total number of moles ” gives 


S 
Seer Cea eee 
x, ee oe eee fe oy ees) (19) 
where X are mole fractions. By definition 
C $ 
a Xt Xx, = 7 (20) 
Jje acet 
Equation (19) is summed over J to give 
C s 
At ye XxX» = zs 
J ee: 4 n (21) 


where 


Substituting for Xi from Eq. (20) into Eq. (21) results in 


S$ 


22 
p+) (y-1)X, =F oe 


Therefore, Eq. (19) may be written 


$s 


oo a “(a Le, 9, H-OYX, Gathered) (23) 


giving © of the S equations necessary for determining the equilibrium compo- 
sition at a specified temperature and pressure. The remaining equations are 
g ee 
= 4 
Xa eX (Leerl, ¢+2,°° 5) (24) 


get 


where Kye are equilibrium constants based on mole fractions. These constants 


are related to Ky, by 


1%-1) 


Kip = Kye CE=c+t, 042, °*'S) (25) 


Kee are calculated from chemical potentials by 


Kpp = exp (- See oe 
s 
where AF; is the change in free energy at standard pressure for the @ th 
reaction. The Afg were calculated by means of polynomial fits to various 
tabulated and derived thermodynamic functions provided by R. E. Duff of 
Los Alamos Scientific Laboratory. Equations (23) and (24) determine the 
equilibrium composition when temperature and pressure are specified; there- 


fore, they can be employed for the reservoir state calculations. 


For computation along an isentropic path, the Brinkley method 
must be modified. The composition and pressure will be computed at various 


temperatures for the entropy at the prescribed reservoir state. Since pressure 


—- —-— mms on Sa | |6 


re to 


is an additional unknown quantity, one more equation is required. A statement 


that entropy is a constant is the desired equation. 


i ~ 2 
tie A toe ie hy En he = tnsze" (27) 


S, is the entropy (per unit mass) calculated at the reservoir state, and 4% 


is the mixture molecular weight which is calculated from 


J = 2 X; Mi 
te 
Sy is the molar entropy of the ith species at standard pressure. Simultaneous 
solution of the nonlinear algebraic equations (23), (24), and (27) gives the 
composition and pressure at various temperatures along an isentrope charac- 


terized by S, . The Newton-Raphson method as employed in the solution is 
outlined in Appendix A. 


Additional quantities are computed from the results of the solution. 
Concentrations in moles per unit mass 7 are found from 
Xe ; 
mes OE ER) (28) 
# 
Density is found from Eq. (4), enthalpy from Eq. (6), velocity from Eq. (2), 
and mass flow per unit aream /A ; from Eq. (3). The state at which mass 
flow per unit area goes through a maximum defines the conditions at the 
throat of the equilibrium nozzle. The local to throat area ratios are then 


found by dividing the maximum m ‘TA 7 by the local m’/A’ values. A Mach 


number is calculated at each step from the expression 


M = u/ Zt) (29) 


(sentropic 
Cre orti an, 


The thermodynamic state is independent of the form of 4(x) because the 


satisfaction of equilibrium conditions at any area is independent of the history 
of the flow. 


The equilibrium air expansions were calculated on the CAL IBM- 
704 computer™ for reservoir temperatures of 4000, 5000, 6000, 7000, and 
8000°K and reservoir pressures of 100, 300, and 1000 atmospheres. The 
calculations assumed a 13 chemical specie system for air. The solutions 
for the gas dynamic properties and gas composition are presented for each 
equilibrium expansion in Tables 1-16. In addition, the air composition is 


shown as a function of area ratio for each case in Figs. 2-17. 


“The Fortran programming for IBM machine computation was done by D. B. Larson, 
Computer Appfications Branch, Systems Research Department, Cornell Aeronautical | 
Laboratory. The IBM-704 program used in the solution of the equilibrium flow 
is availabila at CAL. In the event further equilibrium flow calculations are 
required, details will be provided upon request. | 
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Ill APPROXIMATE SOLUTIONS FOR FINITE-RATE AIR EXPANSIONS 


A. Relaxation Length Griterion for Freezing 


The solution of nonequilibrium: flows, which involve finite reaction 
rates, is rendered difficult by the need for simultaneous solution of inter- 
related chemical rate equations and the gas dynamic equations. Also, at 
the present time, the magnitudes of the required rate constants, and their 
temperature dependence, are still not known with certainty. It is worthwhile, 
therefore, to investigate a simple method for the approximate calculation of 


the final state of a gas undergoing expansion with finite reaction rates. 


The present approach is that of Ref. 3 in which specific non- 
equilibrium airflows in hypersonic nozzles were calculated on the basis of 
a relaxation length. The results of the previous studies“?> have shown 
that for flows involving a single relaxing species, the use of a local relaxation 
length, evaluated for the corresponding equilibrium flow, provides an excel- 
lent estimate of the final frozen state of the gas obtained from exact finite- 
rate solutions. Accordingly, a relaxation length is introduced in the air 
calculations on the basis of an assumption concerning the relative importance 


of the various chemical reactions. 


The energy associated with the normally inert degrees of freedom 
for the air species, as a fraction of the total enthalpy, is shown in Fig. 1. 
It is seen that over the range of reservoir conditions pertinent to the present 
frozen expansion calculations, i.e., from 4000 to 6000°K, oxygen dissociation 
is the most important internal mode as it represents the source of greatest 
energy loss due to freezing. The percentage of the gas enthalpy associated 
with nitrogen dissociation is small compared with that of oxygen. Furthermore, 
since the energy associated with the formation of nitric oxide is small and 
essentially constant over a wide range of temperatures and pressures, the 
potential energy loss due to freezing of the nitric oxide composition is also 
quite small. Consequently, a relaxation length is derived assuming that only 


oxygen dissociation-recombination reactions occur at finite rates, 


11 


The reaction may be written 


c Cé) 
’ é: ; 
OE Onl <0 O2 + Me 
R 
where M, = oxygen atom 
M, = oxygen molecule 
M, = nitrogen atom 
M, = nitrogen molecule 
Me = nitric oxide molecule 


The equilibrium constant in terms of molar concentrations per unit mass 


x) has the assumed form 


ae be Ss 
K. (T ) = ae re x “ (30) 
R (Are 


(o : ae ; 
where # poe reaction rate coefficients. The subscript e denotes the 


local equilibrium value. The rate equation is then 


ae, an, (2) 4,2 () 3 (2 pf &) y»2 
aie eu EY = ve ‘k %, + 2 "bk ss fe %, a af tee X% %, 
(31) 
r2 ¢f) LY Ci) fa 2 os ti) 
-Lp Re A t Zn %, Re On, 4p a Re mM; 


On substituting Eq. (30) and rearranging, Eq. (31) becomes 


? 2 = 
r dh, ’ 2 2 Wa Be &% yi Gore 
u a tev K, k ¢ %, 4 fe a t+ peel Mg 
ha Gees : Ke %, Bs %, ¢73 es %, (32) 


12 


where 


4 Cc) 
apy he 
é kp 


Now the conservation equation for oxygen atoms is 


a aes aie 2%, = CONST. 


therefore h, gis a where 8 = ConsT See 
2 


Defining 


> 
pti 
tt t 
Ra 
NP 
eS 
=| 
hb 


= £ = half throat height (cm) 
Q. tangent of nonal semi-angle 
ae 
vy = Z £0, we (2) 
ES toast 
pr 
and substituting in Eq. (32), noting that K, (Eq. 30) is defined for equilibrium 


values of the 7: , there results 


: 


dx, : 
oa - EA 1 +A, oo $ By at g- & 
uke 9 To ca Yo, é (33) 
BYie 
%, Rs B n a, 7 Io6 
seat} 
\ 2 
where J, = local finite-rate equilibrium value of 25, i.e., the equilibrium 
de q 


value of gj for the actual temperature and density distribution. 


Equation (33) is then of the form introduced by Heir’? 


TMes 5 
(34) 
Cae eee 6 
poe 2 
where 
ta uOx 
ae 5, © om x, ET aa ae 
2 0 M; jose) oe 
Le eas Reiss B ere 


The rate equation, Eq. (34), which incorporates a local relaxation 
length *, is used to derive an approximate freezing criterion for the finite- 
rate flow. The freezing criterion is given ina general form in Ref. 3 and its 


development is repeated here for completeness. 


For flow near equilibrium, J, -%,. << %, and it follows from 


Eq. (34) that ae | ae Joe . For near frozen flow, 7), << J, so that 

4 Yo Re ge. >> %, =, Thus anindication of where significant freezing 

Ax r r oy x 

has occurred may be obtained by setting ae a This equality 
x r 


may be evaluated on the basis of the infinite rate equilibrium flow by using the 


average values of LY, and %, through the freezing region. That is 


x 
at -+ 
A*¥ A NAR 


Ne = he 


(g 


to give 


au.) eile os (3) (= as () 
Ae Jy Pa bE ie he 


where the subscript o refers to infinite-rate equilibrium values. 


Some remarks should be made at this point on the selection of 


"equilibrium" values. As mentioned earlier in this section, the concept of 


14 


a local chemical relaxation length has been found useful for flows in which 

only a single mode of relaxation is involved, such as the dissociation relaxation 
of a diatomic gas, with or without inert diluent. For the more complex case 

of air, the relaxation length criterion has been applied to two simplified 

kinetic models of air in which the only finite-rate reactions considered are 

the oxygen dissociation and recombination processes. All species except 

O and O, were assumed to act only as second or third body colliders for 

the oxygen dissociation- recombination kinetics. The appropriate reaction 

rate constants were obtained from simple collision theory. Vibrational 


equilibrium was assumed for all molecular species. 


Two basic models were considered in the calculations. Inthe 
first, the nitrogen atom and molecule and nitric oxide mass concentrations 
were frozen at their respective reservoir values during the expansion. In 
the second, these mass concentrations were taken as those for infinite- rate 
equilibrium (Sec. II). The flows with oxygen atom-molecule equilibration 


corresponding to these two models were used in the solution of Eq. (36). 


Both models gave values for the frozen degree of oxygen dissoci- 


ation (~p ) which agreed quite closely up to stagnation temperatures of 6000°K. 


For example, in a typical hyperbolic nozzle (4/a = 1.0) for a stagnation 
temperature of 5000°K and a stagnation pressure of 100 atm., Xp = .125 
with N, , N , and NO frozen at their reservoir concentrations and dp =. 149 


with N, , N , and NO in infinite-rate equilibrium. Since these are the 
limiting conditions for the N. ,N , and NO reactions, this comparison 
demonstrates that oxygen dissociation-recombination reactions play the 
major role in air kinetics. The second model has been used in the calcu- 
lations reported herein, i.e., the equilibrium solutions of Sec. II were used 
to evaluate relaxation lengths from Eq. (35). 

In the above calculations the magnitude of ne at 3000°K has been 
taken as 1.16 x 10° em® mole”? sec, from the results of Byron's® inter- 
ferometric studies of oxygen dissociation rates. Then since the temperature 
dependence of the recombination rate coefficient has been assumed as 


RO =e 77 (following Ref. 3) 


I5 


then 
6 
(2) 15 3000 \* cm 
tee, = 1.16x 10 e ——s— 
ie mole™ sec 


The remaining rate constants for the second or third bodies (4;) were 
determined from ordinary collision theory. The form of the rate coefficient 
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is, from collision theory 


ta 


KOE (& ii a 
i ee EE 


whe re K = steric factor 
Dp’ = average diameter of colliding particles 


Aij = reduced mass of colliding particles (cand j) 


_ Ait Aj 
AA 65 
Lg = constant 
o- = symmetry number = 1 unlike particles 
= 2 like particles 
oe = characteristic dissociation temperature for species A 


For oxygen (from Ref. 8), 


(f) ie 
Ea =f, - ois 
Re 2, 


For A, and NO the steric factors and molecular diameters were assumed 


the same as for O, , whence 


Cf) = 
Re, Tenn 


(2) 
Rp PAST, Oy 


16 


a= oom oo Oe Oi eit eee 


then for N, ; 


prs +A, + 2.07 


4 (5) 
2,7 = 45 =O SS 


ForA , the steric factor and atomic diameter were assumed to be the same 
as for 0, so 


T ig 
mee) ee 
ke” ke? fey (2) 


95 


= fe 
1.047 (3547) ~ 36.6 


Relaxation lengths may then be calculated from Eq. (35). 


Equation (36) is solved graphically to locate the point of freezing 
in the nozzle. The slopes pra B are obtained graphically from tangents to 
the curve of oxygen atom mass concentration as a function of area ratio. 
Calculated values of (% /r), and the value of (fh) dA are plotted over 
a range of area ratios. The area ratio at which these two curves cross is the 
location at which significant freezing has occurred. The frozen value of 7% 
that can be expected for finite-rate flow is given by the value of Do 0 which 


corresponds to this area ratio. Figure 18 shows the solution of Eq. (36) for 
the 5000°K, 100 atm. case. 


B. Nozzle Geometry Considered 


The nozzle geometries considered are given by the general expres- 
sion A = +x" These nozzles are symmetric about the throat for even values 
of the exponent 7. ln particular, two nozzle geometries have been considered 


which correspond to values of the exponent equalto land2. For ” = 1, the 


17 


nozzle is a simple two-dimensional wedge-type nozzle with a sharp throat. 
The nozzle parameter, a Lfe ,» which enters the equation for the area 
ratio through the relationship ~ < = , is given by one-half the throat height 
in centimeters divided by the tangent of half the total expansionangle. For 
nm = 2, the nozzle is an axisymmetric hyperbolic nozzle with a smooth throat 
and which becomes a conical nozzle for large values of the area ratio, i.e., 
for nozzle areas not close to the throat. The nozzle geometry parameter is 
given by the throat radius in centimeters divided by the tangent of the asymp- 
totic half angle. 


The nozzle geometry parameter, L/a@ » represents the rapidity of 
the nozzle expansion. For a given throat size, small values of L/a corres- 
pond to large expansion angles and therefore to rapidly expanding nozzles. On 
the other hand, large values of the parameter L/a. represent small expansion 
angles and accordingly correspond to gradually expanding nozzles. Many of 
the nozzles currently used for hypersonic flow research are designed with a 
value of Lfa. of about unity. In the present analyses, values of Lfa of 0.1, 


1, and 10 have been used to bracket the values of current design practice. 


C. Results of Approximate Finite-Rate Solutions 


The results of the relaxation length criterion calculations for flow 
in a hyperbolic axisymmetric and wedge-type nozzle are shown in Figs. 19, 
20, and21. The calculations were performed for stagnation temperatures 
of 4000, 5000, and 6000°K and stagnation pressures from 100 to 1000 


atmospheres, for three values of the nozzle geometry parameter ie. 
(O15. G7 10); 


Figure 19 shows the area ratio for freezing (Ag ) plotted against 
the stagnation pressure for each stagnation temperature. ln all cases freezing 
occurs fairly early in the nozzle (Ar < 25). Further, freezing in the nozzle 
is delayed by an increase in reservoir pressure, an increase in reservoir 


temperature, and by the use of gradually expanding nozzles (larger L/ay, 
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The frozen degree of oxygen dissociation, -¢ (the mass fraction of oxygen 
8 8 f yg 


atoms to total oxygen) is given by 
vo 
Ze Tg fo eZ, Wong 
where the subscript f denotes the value at oxygen freezing. The fraction of 
the stagnation enthalpy represented by frozen chemical energy is given by 


9 je) 
FELT Pr 


Ho Ho 2 


H Eye} 9 
a me SB Flteaah go, 


where 


f 
The characteristic dissociation temperatures (@»o; ) assumed for oxygen, 


nitrogen and nitric oxide are 59,400°K, 113,200°K, and 75,500°K rkanectiels coe 


The frozen degree of dissociation and He / Ho are plotted versus the 
stagnation pressure for each stagnation temperature in Figs. 20 and 21 respec- 
tively. Figure 20 shows that the frozen level of oxygen dissociation is 
appreciably reduced by an increase in stagnation pressure. This indicates 
that high stagnation pressures are required in order to reduce the static 
enthalpy loss through freezing at high levels of dissociation. Figure 21 shows, 
for example, that at a stagnation temperature of 6000°K and a stagnation 
pressure of 100 atm., as much as 20% of the total enthalpy may be frozen 
chemical energy. A wedge nozzle is superior to a hyperbolic nozzle of the 


same L/o. because of a lower frozen degree of dissociation. 
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IV SOLUTIONS FOR EXPANSION AFTER FREEZING 


For the simplified kinetic model of air, the exact finite-rate solution 
is approximated by assuming equilibrium flow up to the freezing point followed 
by a chemically-frozen flow downstream. This approximation was introduced 
by Bray! for the case of a pure diatomic gas. The area ratio for freezing as 
determined by the methods outlined in the last section, represents a unique 
point in the equilibrium solution. The equilibrium gasdynamic properties and 
composition at this area ratio are the initial values for the subsequent chemi- 


cally-frozen expansion. 


In the frozen expansion calculations, vibrational equilibrium is assumed 
with the energy taken as that for a quantized linear oscillator. The total gas 


enthalpy for translational, rotational and vibrational equilibrium, may then be 
written as 


e 
Hos 3T + M4 he (z a ) + My, i oe 
ales @ %ma/Ty 
(37) 
Ow, PN 
+ On, ie “ng + % (7 + SP.) 
+N (7 + ys) + Sao (Pa > Pena : *oe) | 
2 2 
where & = molecular weight undissociated gas (28.85) 
Cy: 
Oy; = ee 


Q 


The characteristic vibrational temperatures ( ey; ) used for oxygen, nitrogen 
and nitric oxide are 2300°K, 3390°K and 2740°K, respectively. 


The method of solution is as follows. Since the gas composition is now 
frozen, all 7; in Eq. (37) are known constants and the gas enthalpy may be 


simply determined for any specified temperature ratio T. The velocity is then 
calculated from Eq. (8). 
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The frozen expansion is isentropic (by virtue of zero reaction rates 14) 
hence 
a 
dH- “F -0 
a (38) 


Then from Eq. (10), in which poms. d, and i J; are now constant, we have 


yx Loy 
- tri = try = ConsT. 
e rE Os ge (39) 


where) y=) J; at freezing 
(5. oe eas ee th )e 


Substituting Eq. (37) for H, the left side of Eq. (39) is evaluated to obtain 


f Sey pay 
de(2 Bh Ln BP an 
ie 2 AP % 
6 
LA 1 Gy, any e % py 
if Oy /T if 6, ta 
et e@ % / e a 
Cy, ‘ 
y die ee ie 7 one 2/74 
7 ( LE an roy a Ov, fT 
ao \ bs Be er 
‘ By {T. 
g 1 ov wie os), 
+ LN ein NO ae oY, 
oy (+ a ne /T 4 . at): mat inp = CONST, 
(40) 
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where 


i 


de 
pie 


Equation (40) relates 7 and T explicitly for flow with frozen dissociation and 


Ga i“ Ne + ONo yp, 
i hae he 


4 


equilibrium vibration. The constant is evaluated initially for equilibrium values 
at the point of freezing. The density is then calculated from Fq. (40) for all 


subsequent | . The pressure is next obtained from Eq. (10). Also the area 
ratio from 


A = (Au ex 


LO th 


and a Mach number is calculated at each step from 


M = uf EE), 
(ad, ise ntropcc 


frozen 


The frozen air expansions were calculated on the CAL Datatron computer* 
and the results for all initial reservoir states and all nozzle configurations 
are tabulated in Tables 17 to 26. The parameters in each table commence 
with their values at freezing, i.e., the equilibrium values at the appropriate 
Ar . The behavior of the flow temperature, pressure, density, velocity and 
Mach number during expansion through the hyperbolic-type nozzles from each 
initial reservoir state is shown in Figs, 22-31. The corresponding equilibrium 
solutions are included to show the calculated location of freezing in the nozzle 
and to provide a basis of comparison for the frozen expansions. Only the 
results for the extreme nozzle shapes (2s 0.1, 10) are shown; data for the 


€=1,0 hyperbolic nozzle flows, as well as for all the wedge-type nozzle flows 
may be extracted from the tables. 


*The programming for Datatron machine computation was done by V.L..Widler, Head, 
Data Handling Section, Hyparsonic Tunnel Dapartmant, Cornall Aaronautical 
Laboratory. 
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Vv FURTHER KINETIC CONSIDERATIONS 


In the application of a relaxation length freezing criterion for air it has 
been assumed that the freezing of the oxygen kinetics, taken here as the gov- 
erning chemical mechanism, implies the simultaneous freezing of all kinetics. 
Actually this assumption has little significance in the present cases since at 
the predicted freezing temperatures, little dissociation energy exists in the N 
and NO, and consequently, only a small amount of enthalpy is involved whether 


these species are assumed to freeze or not. 


However, an analysis similar to that outlined in Sec. IIl may be performed 
for a different chemical species. The appropriate equations would all be of 
the same form as those given previously for oxygen. Now if relaxation length 
considerations were to indicate freezing of this species at area ratios earlier 
than those predicted in Sec. III on the basis of oxygen kinetics, the previous 
assumption of simultaneous freezing of all species would result in a greater 
loss of static enthalpy because of the freezing of oxygen at higher temperatures. 
This in turn could result in significant changes in some of the gas dynamic 
properties calculated in Sec. 1V. lt was considered expedient, therefore, to 
indicate the approximate solution to a specific case as dictated by the kinetics 


of some species other than oxygen. 


The 6000°K, 100 atmosphere reservoir case was chosen as it involved 
the greatest amount of initial oxygen dissociation and is therefore the case 
most sensitive to earlier oxygen freezing. Furthermore, since the equilibrium 
solution exhibited a more rapid decrease in 4, compared to 7, (Fig. 9), the 
analysis was based on the nitrogen recombination kinetics. Citing only the 
results, the application of a relaxation length criterion indicated the freezing 
of nitrogen approximately at the throat for all nozzle configurations. Then, 
consistent with the previous assumption, the nitric oxide and oxygen were 
frozen simultaneously with nitrogen and the resultant frozen expansion is 


indicated by the dashed profiles in Fig. 29. 


lt is seen that the calculated frozen expansions determined Oy the nitrogen 


kinetics may differ appreciably from those determined by oxygen kinetics. The 
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difference is mainly dependent on the location of oxygen freezing. However, 
an argument against the freezing of nitrogen is provided by the nitric oxide 


"shuffle" reactions. The reactions 


N+NO + Ne+0 


Siow 


N+0, —s= WNO+O 


provide two fast two-body processes for equilibration of V- N, and V- NO 
which have not been considered in the nitrogen three-body processes discussed 
above. 
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VI CONCLUSIONS 


A method has been described to estimate the effects of finite chemical 
reaction rates on the one-dimensional expansion of air in hypersonic nozzles. 
Numerical solutions have been computed for reservoir temperatures from 
4000°K to 6000°K and reservoir pressures from 100 to 1000 atmospheres. The 
calculations have considered a range of wedge and axisymmetric nozzle shapes. 
The air was assumed to be at rest and in chemical equilibrium at the reservoir 
state. The flow parameters and gas composition in the infinite-rate equilibrium 
expansion from each reservoir state were calculated, and a relaxation length 
criterion, based on the equilibrium solution, was used to indicate the approximate 
location of freezing for finite~rate airflow in each nozzle configuration. The 
resultant frozen air expansions have been presented in tabular and graphical 


form. 


The results of the analysis have showr. that in all cases, freezing occured 
fairly early inthe nozzle downstream of the throat, at area ratios less than 
about 25 (Fig. 19). Asa result of such early freezing the frozen degree of 
dissociation may be quite large (Fig. 20) representing an appreciable loss in 
the static enthalpy, It is seen, for example (Fig. 21) that in the high tempera- 
ture, low pressure cases, the unavailable enthalpy owing to freezing may be 


as much as 19% of the stagnation enthalpy. 


The following features are evident in the results. The effect of the static 
enthalpy loss through freezing is to reduce the pressure, velocity and tempera- 
ture at a particular nozzle location (area ratio) from the corresponding equilib- 
rium values, and to increase the density and Mach number. While the change in 
temperature and Mach number may be considerable, the density and velocity 


are relatively unaffected by freezing. 


As expected, the effects of freezing are minimized for larger aL /a. values, 
i.e., by the use of gradually expanding nozzles. In practice, however, limita- 
tions will be imposed by boundary layer growth in the nozzle and a compromise 
will be required between maximizing Lia and maintaining a reasonable nozzle 


size. In this respect the wedge nozzle is superior to the hyperbolic nozzle for 
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a given L/o by virtue of a lower frozen degree of dissociation of the flow 
(later freezing). It should be noted that the effect of a change in the nozzle L/a 
parameter also reflects the effect of a change in the kinetic rates, in particular 
Ly ae Since Lio and the recombination rate constant tea. occur in the rate 
equation as a product (in p ), a change in Lfa. is comparable to a similar change 
in ee, the product remaining the same. This factor permits a simple com- 
pensation for discrepancies in existing rate constant data. Also noticeable in 
the results is the reduction in the frozen degree of dissociation with increase 
in reservoir pressure. This indicates that high reservoir pressures are 
required in order to reduce the static enthalpy loss through freezing in the 
operation of hypersonic wind tunnels at high stagnation temperatures. This 
requirement is compatible with aerodynamic requirements for testing at very 
high Mach numbers. 


The present nonequilibrium airflow calculations are meant as approxi- 
mate estimates of the gross features of particular finite-rate nozzle flows. In 
the interests of simplicity several effects have been omitted from the approxi- 
mations used, notably the neglect of the nitrogen and nitric oxide reactions at 
finite rates in the simplified kinetic model of air. However, since the 
assumption of both zero and infinite rates for reactions involving these species 
resulted in relatively small changes in the frozen levels of dissociation 
(Sec. 3.1), the calculated expansions from each reservoir state should repre- 
sent realistic approximations to actual gas behavior in the nozzle. Ultimate 
verification must come from comparisons with exact sulutions to more com- 


plex kinetic models. Such a program is presently underway at CAL. 
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APPENDIX A 


NEWTON-RAPHSON METHOD!?!:!2:13 por 


SOLVING THE EQUILIBRIUM PROBLEM 


The use of the Newton-Raphson method for the iterative solution of these 
simultaneous nonlinear algebraic equations was originally suggested by 
Brinkley” for the determination of equilibrium composition at a given temper- 
ature and pressure. Brinkley's method is directly applicable to the computation 
of the stagnation state; however, certain modifications are introduced for 
computing composition and pressure at a given entropy for various temper- 
atures along the expansion. Only the modified form is presented here because 
the essential features of the Rrinkley method are contained therein. Equations 


(23) and (27) are expressed as 


i =O Cy ci ee Seon a Pao oh) (A-1) 
where 
S ‘ 
& rye sya | 2, (% a4) i = 4; | X¢ oA (A-2) 
alot = 
(Raya 
he De k: (% - My So - tnX; ) - trp’ (A-3) 


For given J, and rage , values of OS and p- (the (C+? )st independent variable) 
must be guessed such that Equations (A-1) are satisfied. xX; are related to 
x; by Equation (24). 


ies yy. 
pF acm Goad 5 ae) 
Iai Sime Cre fe? es) eat aa) 
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The Newton-Raphson method provides a means of improving successive 
guesses by utilizing a Taylor expansion of Fy about xy and p to first order 
terms. This leads toa set of fractional corrections hg” (kh <?, PRES E a He at | 
applicable to the xy and p values obtained at the rth iteration. The improved 


(+1? )st values of the independent variables are found from 


(ret) Cr) (r) , 
LSS PEGE Tee) Chet, -+rc) (A-4) 
and 
1st) ate? (r) 
a CRY. od 


The set of hp is computed from a set of linear equations expressed in matrix 
notation by 


lag’ | [42] - | 


hy, and Fare column vectors 


) 
Fi 2] Cit, 2, ee Cis re (A-6) 
CRE 2, eee at) 


£ 
Aik = Xj Shp ee Ye | cs Gi C) a 1) X¢ (A-7) 
3; =O for jx Cf atye2. oe) 
5; @ efor 4 =k Chk=i,2,--+¢) 


Ajsens * (4-1) |r; 47 O%y-1)) Xe fet a,c) (A-8) 


vacel 
Aest,4 ~[1-(S4 “Ak So oa dn X,) | Xx + 71 
é A P A-9 
i, Veg [7- CS, ~My So - nity) | Xz i (A-9) 
? *1,2,-"6) 


Acstenm 1 + yi (Me -1)|1-(G ~My 55 - On Xy)) Xe (A-10) 
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The (7 +?/)st iteration is started by calculating hp from Equation (A-6) using 
the formulas in Equations (24), (A~2), (A-3), (A-7), (A-8), (A-9), and (A-10), 
Improved values of Xj and p. are found from Equations (A-4) and (A-5), and 


é 0 : * , t 
improved values ofX4 are found from Equation (24) using X ; (r+?) ana p aes I 


| 
The cycle begins again with the computation of [A 1 | from Equation (A-6). 
When all hy become equal to or smaller than a prescribed maximum fractional 


error €, , the computation is terminated. 
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Figure 21 (b) EFFECT OF STAGNATION PRESSURE AND NOZZLE GEOMETRY ON THE 
FROZEN CHEMICAL ENERGY FOR A SIMPLIFIED AIR MODEL 
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Figure 21 (c) EFFECT OF STAGNATION PRESSURE AND NOZZLE GEOMETRY ON THE 
FROZEN CHEMICAL ENERGY FOR A SIMPLIFIED AIR MODEL 
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